N AS A/CR-2002-2 11924 
ICASE Report No. 2002-34 



Linear Parameter Varying Control for Actuator 
Failure 

Jong-Yeob Shin 
ICASE, Hampton, Virginia 

N. Eva Wu 

Binghamton State University of New York, Binghamton . New York 
Christine Belcastro 

NASA Langley Research Center, Hampton, Virginia 



October 2002 



The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated 
to the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key 
part in helping NASA maintain this 
important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the lead center for 
NASA’s scientific and technical information. 
The NASA STI Program Office provides 
access to the NASA STI Database, the 
largest collection of aeronautical and space 
science STI in the world. The Program Office 
is also NASA’s institutional mechanism for 
disseminating the results of its research and 
development activities. These results are 
published by NASA in the NASA STI Report 
Series, which includes the following report 
types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results 
of NASA programs and include extensive 
data or theoretical analysis. Includes 
compilations of significant scientific and 
technical data and information deemed 

to be of continuing reference value. NASA’s 
counterpart of peer-reviewed formal 
professional papers, but having less 
stringent limitations on manuscript 
length and extent of graphic 
presentations. 

• TECHNICAL MEMORANDUM. 

Scientific and technical findings that are 
preliminary or of specialized interest, 
e.g., quick release reports, working 
papers, and bibliographies that contain 
minimal annotation. Does not contain 
extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATIONS. 
Collected papers from scientific and 
technical conferences, symposia, 
seminars, or other meetings sponsored or 
cosponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to 
NASA’s mission. 

Specialized services that complement the 
STI Program Office’s diverse offerings include 
creating custom thesauri, building customized 
data bases, organizing and publishing 
research results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home 
Page at http://www.stLnasa.gov 

• Email your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA STI 
Help Desk at (30 1 ) 62 1-01 34 


• Telephone the NASA STI Help Desk at 
(301)621-0390 

• Write to: 

NASA STI Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076-1320 



N AS A/CR-2002-2 11924 
ICASE Report No. 2002-34 



Linear Parameter Varying Control for Actuator 
Failure 

Jong-Yeob Shin 
ICASE, Hampton, Virginia 

N. Eva Wu 

Binghamton State University of New York, Binghamton, New York 
Christine Belcastro 

NASA Langley Research Center, Hampton, Virginia 
ICASE 

NASA Langley Research Center 
Hampton, Virginia 

Operated by Universities Space Research Association 



Prepared for Langley Research Center 
under Contract NAS 1 -97046 


October 2002 



Available from the following: 


NASA Center for AeroSpace Information (CASI) 
7121 Standard Drive 
Hanover, MD 21076-1320 
(301)621-0390 


National Technical Information Service (NTIS) 
5285 Port Royal Road 
Springfield, VA 22 161-2 171 
(703) 487-4650 



LINEAR PARAMETER VARYING CONTROL FOR ACTUATOR FAILURE 

JONG-YEOB SHIN*, N. KVA WUp AND CHRISTINE BELCASTRO* 


Abstract. A robust linear parameter varying (LPV) control synthesis is carried out for an HiMAT 
vehicle subject to loss of control effectiveness. The scheduling parametc r is selected to be a function of the 
estimates of the control effectiveness factors. The estimates are provided on-line by a two-stage Kalman 
estimator. The inherent conservatism of the LPV design is reduced through the use of a scaling factor on 
the uncertainty block that represents the estimation errors of the effectiveness factors. Simulations of the 
controlled system with the on-line estimator show that a superior fault-tolerance can be achieved. 

Key words, fault tolerant control system, fault parameter estimation, reconfigurable controller 

Subject classification. Guidance and Control 

1. Introduction. One of control schemes for a nonlinear system is a gain-scheduled linear parameter 
varying control technique [13, 1. 7, 16]. This approach is particularly appealing in that a nonlinear plant is 
treated as a linear parameter varying (LPV) system whose state-space matrices are functions of a schedul- 
ing parameter vector. This allows linear control techniques to be applied to a nonlinear system. Several 
researches on an LPV synthesis methodology allow the design of the global control law for an LP\ system 
over a parameter set which is bounded and measurable [13, 1, 7. 16 . An LPV controller synthesis is formu- 
lated into a linear matrix inequality (LMI) optimization problem. Then are LPV control synthesis methods 
according to a functional form of an LPV system. The polytopic LPV control synthesis method [2] is used 
for an LPV system which is a polytopic function of a scheduling parameter vector. The affine LPV control 
synthesis method [1] is applied to an affine LPV system, whose LMI constraints are evaluated at only vertex 
points of an LPV system. The grid LPV control synthesis method [7, 16] is for an LPV system which is a 
bounded function of a scheduling parameter vector. In the method, LMI constraints are evaluated at grid 
points over parameter spaces. These methods can be converted to each other by increasing conservatism 
to describe an LPV' system. The grid LPV synthesis method has been successfully applied to synthesis 
controllers for the pitch-axis missile autopilots[17, 12], F-14 aircraft lateral-directional axis during powered 
approach [6, 4], turbofan engines [15, 5] and F-16 aircraft [14]. Scheduling parameters of these applications 
are physical parameters such as angle of attack, mach number, velocity, dynamic pressure, etc. Schedul- 
ing parameters in LPV control synthesis are required to be measurable and the variations of scheduling 
parameters should be in a bounded set. 

In this paper, actuator failures are modeled as an LPV system as functions of actuator effectiveness 
parameters [9]. These parameters are estimated as biases using an augmented Kalman filter. A set of 
covariance-dependent forgetting factors is introduced into the filtering algorithm. As a result, the change in 
the actuator effectiveness is accentuated to help achieve a more accurate estimate more rapidly. The 
bounds on parameter estimation errors are assessed through simulations, which are then used as bounds 
of real parameter uncertainty in the construction of a robust LPV control law. Actuator faults can be 
parameterized as estimated fault effectiveness parameters. Thus, it is possible to formulate a fault tolerance 
control design problem as an LPV control synthesis problem based on estimated faults parameters. 

Fault estimation errors and modeling uncertainties are represented by an uncertainty block in the con- 
struction of a robust LPV control law. The structure of an uncertainty block is not included in an conventional 
LPV control synthesis methodology [7, 16]. A scaling factor on a uncertainty block can reduce conservatism 
of the LPV synthesis [1], In Ref.fl], it is formulated into a single optimization problem to find a scaling 
factor and a control law to achieve a certain level of performance. The optimization problem is not a convex 
problem, which has unknown positive matrices X and Y related with a control law and a scaling factor 5 
related with the uncertainty block structure. The problem is solved by m iteration method of fixing X and 
Y or a scaling factor S. In this paper, the problem is formulated into two LMI optimizations: one is to design 

’ICASE, MS 1320, NASA Langley Research Center, Hampton, VA 23681. This research was supported by the National 
Aeronautics and Space Administration under NASA Contract No. NAS1-97046 while the first author was in residence at. 
ICASE. NASA Langley Research Center, Hampton, VA 23681. 

^ Professor, Department of Electrical Engineering, Binghamton State L Diversity ol New York, Binghamton, New \ork, 13902. 

I Researcher, Guidance Control Branch, NASA Langley Research Center, Hampton, VA 23681. 
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an LPY controller with given a scaling factor and the other is to determine a scaling factor 5 with a given 
control law. The problem is solved by an iteration method of fixing a control law instead of fixing variables 
A r and Y. This helps to find a scaling factor 5 for minimizing an induced-^ norm of the closed-loop system. 

This paper contains the following sections. In Section 2, an LPY synthesis methodology used in this 
paper is summarized. In Section 3, fault parameter estimation methods are presented. In Section 4. an LPY 
controller for an aerospace vehicle is designed by using an LPY synthesis control methodology with a scaling 
factor. In Section 5, the simulation results of the closed-loop system are presented and this paper concludes 
with a brief summary in Section 6. 

2. LPV Synthesis. 


2.1. Problem Statements. In this section, a control synthesis problem is defined, based on an es- 
timated parameter vector p G 7Z Ua such as actuator effectiveness [9]. Actuator effectiveness parameters 
represent actuator failure cases (actuator damage). Suppose actuator effectiveness parameter can be esti- 
mated using the estimation methods presented in Section 3. A system dynamics can be represented by an 
LPY system according to an estimated scheduling parameter vector and estimation error bounds. 

An LPY system can be represented as functions of an estimated scheduling parameter vector p with an 
uncertainty block A which captures parameter estimation errors and unmodeled dynamic's. An LPY system 
can be written as: 


X 


' A(p) B A (p) B p (p) B u (p) 


X 



C e (p) D AA (p) D Ap {p) Da v (p) 


d& 

e p 


Cp{p) Bp^(p) Dpp(p) Dp U (p) 


d p 

. y . 


Cy(p) Dy A {p) Dyp(p) 0 


u 


( 2 . 1 ) 


where x G 1Z n * . € 7Z n * , e p G 7 Z n *p , y G 7Z Uy , G 7 Z UA , d p G 7 Z nd p , and u G 7l 1lu . All of the state-space 

matrices are of appropriate dimensions. 

An uncertainty block set A is defined as: 


A — {A — diag(S] I \ , • • • , S n I ni A n+ i , • * • , A n + m ) : Si G R, A* G R txl ,cr( A) < /3}, 
where 3 is normalized to 1 without loss of generality. There exists a scaling factor set S such that 

5 = {5 : 5 > 0, SA = AS, S G 7 l nA * n *}. 

The input /output scaling matrices and J 1/ 2 are defined as 


( 2 . 2 ) 


(2-3) 
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^ 'J^(nA + n dp )x (n A +n rfp ) j\/2 _ 
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^(«A + n. Cp )x(n A +n, p )^ (2.4) 


where 5 I//2 G S and 5 G 5. 

The induced £2~nonn °f a parameter dependent stable LPV system is defined as 


| G || 2<-2 


e 2 

SUP TTTj-p 
Vp, d<EC 2 , d ? 0 ||«||2 


for zero initial conditions x(0) = 0. 

Suppose there exists an LPV controller K(p) which stabilizes an LPV system P{p). The control synthesis 
problem is: 


min \\J i / 2 (S)F l (P(p),K(p))L-^' 2 (S)\\^ 2 , (2.5) 

K(p), S€S 

where Fi(P(p),K(p)) means a lower linear fraction transformation (LFT). The optimization problem of 
equation (2.5) is not convex in K{p) and S . The problem is similar to a D-K iteration. In this paper, we 
approach to the problem in a similar manner of solving the D-K iteration problem[3]. 
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2 . 2 . Control Synthesis Methodology. In this section, a procedure of solving the problem is pre- 
sented. There is an LPY control synthesis methodology in Ref. [ 1 ] with a scaling factor. In Ref. [1], an LMI 
optimization can be formulated with unknown matrices A' > 0 . I > 0 , and scaling factor matrices J 1 - 


and L~'/ 2 . However, an LMI optimization problem in Ref. [ 1 ] has an equality constraint. In this paper, to 
avoid an equality constraint, the augmented LPY system with the scaling matrices ./ 1 and L is used 


to design an LPY controller. 

Suppose a scaling factor S is given. The augmented I.PY open-loop system with the scaling matrices 
and can be written as: 



Bi{p) 

Du(p) 

Dnip) 


B 2 (p) 1 x T A{p) Bi(p)L l / 2 

DM d = J'^Ci(p) jWDnip'L-W 
0 J |uj L W) D 21 (p)L -' / ‘ 2 


Blip) 

J l/2 D M 

0 


( 2 . 6 ) 


where c = [c£ c 7 ] r and d - [d' A dj] r . With assumption that D v ,(p) and D-<\{p) are full column and 
row rank for all p. respectively, an LPV control synthesis methodology in Ref. [ 16 , 7 ] can be used in this 
paper. For the sake of completeness, a brief summary of the LPY control synthesis methodology in Ref. [ 16 ] 
is presented in this section. 

Using Q-R decompositions [ 7 ] of matrices D ] 2 (p) and A>i(£), the augmented LPV system is rewritten 


A(p) Bn(i>) B l2 (p) B 2 (p) 

C n (p) D\\\\ ip) £>1112 (p) 0 

Cviip) D\vn(p) B \\22 {p) In,., 

Clip) 0 In d , 0 

where d, £ 7 l Tld ' , d 2 £ II nd * , e { £ R n *\ , and v 2 £ R 71 ** . Matrices Dm,, D [n 2 . £>1121, £>1122, C11 , C v > are 
of appropriate dimensions. 

There exists an LPV controller A (p) which leads to the induced-^ norm of the closed-loop system being 
less than constant 7 . The LPV controller K(p) can be constructed from the solution matrices, X(p) £ R nx ” 
and Y(p) £ R nXn which are calculated by solving the following LMI optimization. 

min 7, (2.8) 

X(/5),V(p) " 

subject to 

X(p)A r (p) + A(p)X{p) + Y.T=i(^)-Bi(p)Bl{p) X(p)CUp) 7 ~ l B{p) ' 

Cn(p)Xip) -In., ~r l D [n .(p) < 0 . ( 2 . 9 ) 

7 ~'B r (p) A-'DJuXp) -I nd . 




Aip)Y(p) + Yip)A T ip) + Z',li “ C Up)Ci{p) Y{p)Bi , (p) 7 "‘C 7 ip) ~ 

B'l\{p)Y{p) - 7 _1 ^n.i(p) <0- (2.10) 

7 - x C{p) 7 -, Cn. ip) -In. . 


'X(p) 7 _1 /„" 

;r l In Yip) _ 


> 0, 


( 2 . 11 ) 


A' ip) > 0, r(p)>o, 

when' 

Aip) = A(p) - Biip)CM, Dip) = Blip) -Biip)D U i.ip), (2-12) 

Aifi) = Aip) - BMCiip), C(p) = Ciip)-D l i- 2 (p)Ci(p). ( 2 . 13 ) 
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The definitions of matrices D \ ]*?., i?n.i, i?ni. and Z?n.2 are taken from Ref. [16]. Also, the realization of 
an LPY controller from the solution matrices A" and Y are taken from Ref. [16]. 

The benefit, of the LPV synthesis methodology is that there is no limitation of an affine functional form 
of LPV system state-space matrices. Since the LMI constraints in equations (2.9)-(2.10) are evaluated at 
grid points over all scheduling parameter spaces, an LPV system should be just a function of a. scheduling 
parameter. In this paper, we consider actuators are failed one at a time. In the failure case, the system 
variations due to actuator failures can not be represented by an affine function of an actuator failure pa- 
rameter vector. The disadvantage of the LPV control synthesis methodology is that robust stability over all 
parameter spaces is not guaranteed unless choosing appropriate number of grid [joints. 

Suppose there exists an designed LPV controller K{p) which stabilizes the augmented LPV system. The 
closed-loop LPV system with a given controller is written as: 


id 


Ael(p) Bi nl (p)S- 1 ' 2 Ihjp) ' 


X c l 

ca 

— 

S'^CiJp) S l ' 2 DnAp)S~ 1/2 S'^DnJp) 


dA 



C->Jp) D 2lrl (p)S-W D T2cl (p ) 


(I p ~ 


where x c \ = [x 7 x^] T . 

Applying the Kalman- Yakubovich- Popov (KYP) Lemma [9], the LMI optimization is formulated to find 
a scaling factor S. There exists an scaling factor S € S which leads to the induced-£ 2 norm of the closed-loop 
system being less than 7*. The scaling factor S can be determined solving the following LMI optimization. 


mm 

r>o,ses 


7 s 5 


(2.15) 


'M 11 

PB 1 

+7.“ 'Cl,Dn cl +Cj cl SDn rl 

PB ic i +7s l Cl ct D‘22 rl + C[ ct SDi2 cl 


★ 

★ 

77 

l Dl rt D 21e ,+DT lfl SDu cl -S 
★ 

i'7 1 dt 1c1 d 22c{ +dT 1c1 sd 12c1 

7s 1 P22 c ,P'X2 ci + Dy2 cl SDi'2 rl ~I 

< 0 


( 2 . 16 ) 


where 


Mu = A r P + PA + P + + Cl,SC irr 

In the LMI constraint, ★ denotes a symmetric component. 

The iteration procedure to solve the problem in equation (2.5) is follows: 

L Design an LPV controller K(p) for a system from the LMI optimization in equation (2.8) with fixed 
S. At the first iteration, 5 is assumed as I. 

2. Solve the LMI optimization problem in equation (2.15) over P(p) and S based on the closed-loop 
system with the designed LPV controller K{p). 

3. Generate an augmented LPV system with the scaling factor 5 


Gh-i(p) = 




■5- 1 / 2 O' 

0 I ’ 


(2.17) 


where Gj(p) is an LPV model at the i th iteration. 

4. Iterate over step 1 to 3 until convergence or terminate iteration based on satisfaction with a designed 
LPV controller. 

The iteration method can not guarantee finding global solutions of K and S since the problem in 
equation (2.5) is not convex in K and S. Also, there is no guarantee of convergence in the iteration process. 
In the LPV synthesis methodology, the matrix P(p) is related with the solution matrices X and Y of 
equations (2.9)-(2.10). When a designed controller is fixed to calculate the scaling factor S, the matrix 
P(p) can be calculated from the solution matrices X and Y and fixed in the LMI optimization [16]. In this 
paper, the matrix P(p) is also set as an unknown matrix in the LMI optimization in equation (2.15) to relax 
the constraints of fixing the LPV controller K{p). Thus, there are two LMI optimizations in the iteration 
process. 
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3. Parameter Estimation. This section briefly describes the formulation of a real parameter estima- 
tion problem, which, when specialized to the actuator effectiveness estimation, transforms a fault (loss of 
actuator effectiveness) tolerant control problem to a robust LPY control problem. The development of this 
section follows that in Ref. [9]. 

The estimator is based on a linear discrete design model of the form: 

' u ). ' 

J * +1 = A ,l k x k + [bf}y l k ••• fell*"] 

. «*" . 


+B k u k + w k (3-1) 

= A d k x k + Efa k + B*u k + (3.2) 

7*+i = 7 * + «>*. ( 3 - 3 ) 

y k = C k x k . + Vk, ( 3 - 4 ) 


where x k G i? n *, G R n " , u k G R Hu and y k G R 7ly are the state, bias, input, and output variables, 
respectively. The discrete model can be obtained from a continuous model via, for example, the Euler s 
rule with a sampling period T s . which preserves the functional dependence of the k 'B matrix on 7'* The 
bias vector 7 with component — 1 < 7 ^. < 0 relates to a actuator failure parameter. It is obvious that 
Ei = B<1 x diag{u [ k , ■ • .u£ u }. w k , uP k and v k denote the white noise sequences of uncorrelated Gaussian 
random vectors with zero means and covariance matrices Q k , Q k and R k , respectively. 

The minimum variance solution is obtained by a direct application of the two-stage Kalman filter al- 
gorithm of Keller and Darouach[ 8 ], with constant coefficient matrices in Ref. [ 8 ] replaced by time- varying 
matrices. The filter is decoupled into four sets of equations. They given as follows. 

- Optimal bias estimator 

7*+i|* “ 7*|* > 

n\.i k =p ii\k+Qi 

7*+i|*+i =7t+i|* + A^' + ,(r*+i — H k fn*-7/r|*)> 

^ V A + l = P k+l\k H k + \\k 

(Hk+i\kPk + \\k H Z+i\k + ^*+ 1 ) '• 

Bias-free state estimator 


**+i|* = ^*^*|* + B d k u k + w k%\k - ^ *+ 1 1*7 *|* 7 

PZ+m = 4^i*(4) T + OS + 

- t *+ 1 1 * R k + 1 1 k ^ *+ 1 1 * » 

^*+l|*+i = ^fc+l|* d" R k +i (y*+i ~ j £*4 1 j* )? 

1 = 1 1 >, ( ^ I ^ 1 1 ^ 

+fi*+l} 1 1 

P k +l\ k +] = (I - Rk+ l^k+l )T*+1|*’ 


where the filter residual and its covariance are given as 


^*+1 — Vk + 1 — 

s k+l = Ct +l PZ +Uk (C d k+l ) r + R k+l 



Coupling equations 

W k = .4*1*1* + E d k , 

Vk + w = w k P2 lb (r^ b r', 
Hk+i\k = C d k+i \ k+llk , 

U+l|*+l = U+l|fr - K k +lHk+l\k- 

And finally the compensated state and error covariance estimates 


%k + 1 \k + \ ~ + ^ 4 +l|A'+l 7 *+ 1|*+1 7 

flfc+l|*+] = Pk+]\k+l + ^Wl|fc+l^* + i|jfc + i^Jb+l|it+l- 

A further measure is taken to modify the above filtering algorithm so that the estimates become more 
responsive to abrupt changes in the control effectiveness factors. 

A well known technique for estimating time- varying parameters is the use of forgetting factors. The 
basic idea is to enable a recursive algorithm to discount the past information so that the filter is more apt to 
recognize the changes in the system. Since the time update of the bias estimate governed by 7 A _ j _ 1 | A — 7 * 1 * is 
the dominant opposing force to acknowledge the abrupt changes in the biases, forgetting factors introduced 
into the time propagation equation P k+l \ k = ^k\k + Ql ^ ie ^ias covariance is likely to function most 
effectively. 

Assume that covariance P k \ k adequately describes the bias estimation error along both temporal and 
spaeial directions under the normal system operation condition. Then this covariance provides a basis for 
the selection of forgetting factors. The bias estimates should be prevented from being impetuous, as well 
as from being indifferent to the changes shown in the measurements. A technique suggested in Ref. [ 11 ] 
amounts to select forgetting factors that would force the adjusted covariance in = P^ k + QJ to stay 

within some prescribed bounds 

G’nmi-f ^ ^A’+ 1 |A- — ^"niax-fj (3.o) 

where rr m i ri and a max are positive constants with 0 < cr m in < ^max < 00 , and I is the identity matrix. Let 
the dyadic expansion of P^ k be given by 


06) 

1=1 

where a[| A , are the eigenvalues of P^ k with > ... > a^, and c|, ..., e k u are the corresponding 
eigenvectors with ||c A || = ... = ||e A u || = 1. Equation (3.6) can then be expressed as 


pi 

r *+i|* 


= 0 < Ai < 1 . 


i=] 


Following the argument in Ref. [11], the forgetting factor A a can be chosen as a decreasing function of 
the amount of information received in the direction e A . Since eigenvalue o A j A of P^ k is a measure of the 
uncertainty in the direction of e \ , a choice of forgetting factor A a based on the above constraints can be 


AL = 


A o 1 

OjbjA: I Omiu + 


L * -- ; . ‘ .n u a Q 

31 max * 1 * 


*-*k\k ^ ^max; 
j ' a k\k — a niax- 


(3.7) 


The estimation algorithm discussed in this section will be seen to have been applied successfully to a HiMAT 
vehicle. 


4. Example. In this section, the iteration process described in Section 2.2 is applied to control a 
HiMAT vehicle for actuator failure cases. Recall it is assumed that actuators are failed one at a time. Thus, 
the control reconfigurability of the HiMAT vehicle never goes to zero [10]. The system variations due to 
actuator failures can be modeled as an LPV system, a function of an estimated scheduling parameter. 
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4.1. Linear Parameter Varying Model of HiMAT. The model of the HiMAT vehicle taken from 
the //-synthesis Toolbox [3] has two inputs: elevons S< and canards ri r : two outputs: angle of attack n in 
radians and pitch angle 6 in radians; and four states: velocity V in ft /sec, angle of attack, pitch rate q in 
rad/sec, and pitch angle. The open-loop model is 


X 


A D 
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'<5; 

y_ 


P D . 


u 

, u = 
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(4.1) 


where 


—0.0226 
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-0.414 

0 
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-2.63 

0 
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22.4 

0 

0 

1 

0 


0 
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0 

1 


. D — 0*2 x 2 ■ 


(4.2) 


A failure parameter vector t — [ti t 2 ]^ is determined by tlie two actuator effectiveness parameters T\ 
and 7*2 of elevon and canard actuators, respectively. Assume that the failure parameters linearly enter in the 
model. The state-space model of the HiMAT vehicle is written as 


x — Ax H- B(r)u , y = Cx , 


where A and C are constant matrices and B{r) = [b\T\ b 2 r 2 ]. Tlie vectors b\ and by are the columns of 

B. The actuator failure parameters can be estimated using the estimation method described in Section 3. 
However, there is estimation error S T = [S Ty 8 r . 2 ] T . The actuator failure parameter vector r is written as 


r = r + S T , 

where f is an estimated value. The estimation error bound is assumed as 8J 8 r < 0.05 for each actuator 
failure case. The matrix B(r) is rewritten as: 


B(t) = B 


T\ 4- 0.05<5i 
0 


0 

T'2 + 0 . 05 ( 5*2 


(4.3) 


where the real uncertainty parameters 8] and S 2 vary from -1 to 1. respectively. 
The LPY model of HiMAT is 


x — A x -f B 


T\ 

0 


0 

T2 


u A B 


0.05 

0 


0 

0.05 




(4.4) 


z = u, y = Cx, 


(4.5) 


w = A 2 , 




(4.6) 


The LPV model in equation (4.4) is a linear function of a parameter vector r. However, the parameter 
vector r can not be chosen as a scheduling parameter since we consider that the actuators are failed one at 
a time. Thus, r\ and r 2 can not be zero simultaneously. To describe the failure cases, a synthetic scheduling 
parameter p is introduced as 0 < p < 2. 

0 < p < 1 : 0 < fi < 1, t 2 = 1 

p = 1 : f\ = 1, f 2 — \ (4.7) 

1 < P £ 2 : T\ = 1 , 0<r 2 <l 

Note that the LPV model of the HiMAT vehicle is not an affine functic n of a scheduling parameter p. 
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Fie. 4.1. Interconnection structure for the model matching problem. 


4.2. LPV Controller Design. The control objective is to track a pitch angle command at actuator 
failure cases. A designed LPV controller should robustly stabilize the HiMAT vehicle over the failure param- 
eter variations. The controller synthesis problem is formulated as a model matching problem in Figure 4.1. 


The ideal response model Tj of pitch angle is taken from the example in the //-synthesis Toolbox [3]. 
The performance weighting function W p and unmodeled dynamics W n are also taken from the example in 
the //-synthesis Toolbox [3]. The sensor noise is modeled as white noise with 0.6° amplitude for angle of 
attack and pitch angle measurements. The weighting functions in Figure 4.1 are 


Tj 

Wp 


s/0.8 + 1 5 
_ . in s/50 + 1 
" 4U s/0.05+1’ 

- n 9 s / 5+1 t 

~ s/iooo + r 2 * 2 ’ 

= 0.01/2x2. 


The control synthesis problem of the HiMAT vehicle is formulated to minimize the induced -£2 norm of the 
augmented LPV system with the weighting functions. 

To solve the control synthesis LMI optimization problem in equation (2.8). basis functions for X(p) and 
Y(p) are required since A" and Y are assumed as functions of p. 

X(p) = J^Mp)X u X, € TZ nxn , Y(p) = ^ i g J (p)Y j , 1' € 7 2"*". (4.8) 

* J 

where basis functions fj(p) and gj(p) are given before solving the LMI optimization in equation (2.8) over 
A^ and Yj. There is no analytic method to choose optimal basis functions for A' and Y in general. The 
functions fi{p) and gj(p) are related with sensitivity of unknown matrices Xj and TJ, respectively. In this 
paper, the basis function set is defined as {1. 1/p, p} for X and Y to help the LMI optimization for total 
failure cases (fi = 0. f 2 = 0). Note that it is not necessary to define that gj{p) is equal to / 2 (p). Since X and 
Y are functions of p, the parameter rate bound is required to solve the LMI optimization in equation (2.8). 
Recall that the scheduling parameter is an actuator failure parameter. Thus, for example, the scheduling 
parameter can suddenly vary from 1 (no failure case) to 2 (total canard failure case). In this paper, the 
parameter rate bound is assumed as \p\ < 100 to capture sudden variations of the scheduling parameter. 

To make the LMI optimization computationally tractable, the LMI constraints are evaluated at the 
following grid points: 


p € {p|0.01, 0.1, 0.2, , 1.9, 2}. 
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Tahi.k 1.1 

7 values in the LMI optimization in equation 2.8) 


Iteration 7 S 


1 

1.23 

diag([l,l,lJ 

L]) 



2 

0.71 

diag( [0.497, 

0.168, 

1.18G. 

1.277]) 

3 

0.60 

diag([1.404, 

1.289. 

1.801. 

1.453]) 

4 

0.54 

diag( [1.580, 

2.193. 

2.100. 

1.187]) 

5 

0.85 

diag( [2.007, 

1.430. 

2.506. 

1.849]) 


Also, the same grid points are used to solve the LMI optimization in equation (2.15). To solve the LMI 
optimization in equation (2.15), the basis function set for P is required Since the matrix P is related with 
A r and Y, the basis function set for P is chosen as {1, \/pJ)}. 

In this paper, the scaling factor S is assumed as constant over all scheduling parameter variations. The 7 
values and the scaling factor S for each iteration are written in Table 4.2. The scaling factor 5 is associated 
with the uncertainty block A which is 


A = diag([8\, <*>2, W)- 


Recall that the real uncertainty parameters rfi and 82 are associated with elevon and canard actuator failure' 
parameters, respectively. The multiplicative uncertainty parameters 8 f ^ and 8 can are also associated with 
elevon and canard control channels in Figure 4.1. 

The iteration process is stopped at the o'* iteration since the 7 value at the iteration is greater than 
the previous iteration. Recall that the iteration process is not guarani eed to be converged. However, the 
performance index 7 in the LMI optimization of equation (2.8) is significantly reduced from 1.23 to 0.54 by 
using the scaling factor 5. In the remainder of this paper, the LPY controller for the HiMAT vehicle denotes 
the designed LPV controller at the 4 th iteration. 



time (sec) 


Fig. 5 . 1 . The LPV controller simulations with and without actuator failures. 
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5. Simulations. In this section, the designed LPV controller is applied to control the HiMAT vehicle 
for pre-defined failure scenarios. One of the failure scenarios is that a total canard failure occurs 1 and 10 
seconds and a total elevon failure occurs between 20 and 40 seconds. The scheduling parameter corresponding 
to the failure 1 scenario is: 



( 0.01, 

1 < t < 10 sec, 


p = < 

1.00. 

0 < t < 1 sec, 10 < t < 20 sec 

(5.1) 


[ 2.00, 

20 < t < 40 sec. 



For the purpose of comparison, the LPV controller for the HiMAT vehicle is simulated both with and without 
actuator failures. The simulation results are shown in Figure 1 5.1. The pitch angle commands are given as 
10° at 1 sec, 0° at 10 sec, and 10° at 20 sec, sequentially. 

In this simulation, the scheduling parameter corresponding to the true failure parameters, as shown in 
the bottom of Figure 5.1 is fed into the LPV controller. It is observed that the LPV controller achieves the 
desired goal of tracking pitch commands in the presence of actuator failures. It can be seen from the second 
and the third plots in Figure 5.1 that the LPV 7 controller always relies on the healthy actuator to track the 
pitch commands abandons the failed actuator. For example, the LPV controller keeps the elevon actuator 
signals close to zero at the elevon actuator failure case. 

The LPV controller is also simulated for the same faulty system as described in equation (5.1) but with 
bounded real parameter perturbations: A p = diag([fi \ , ||A P || < 1. Two examples of simulations with 

perturbations are shown in Figure 5.2. “pertl” and “pert2” in Figure 5.2 denote the cases [S\ S 2 ] = [1 1] and 





KlG. 5.2. Time responses of pitch angle with LPV and H ^ controllers. 

[Ji A>] = [— 1 — 1], respectively. The simulation results show that the LPV 7 controller can robustly stabilize 
the perturbed system and achieve the desired performance level of tracking the pitch angle commands. 

For the purpose of comparison, a fixed H Xl controller is designed at p — 1 (without failures) using 
//-synthesis Toolbox [3] with the same weighting functions described in Section 4. The closed-loop response 
with this H ^ controller is also simulated for the actuator failure scenario of equation (5.1), and shown in the 
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bottom plots of Figure 5.2. It can bo seen that the controller can achieve the desired performance level 
when canard fails(l < t < 10 sec). The controller cannot, however, achieve' the desired performance level 
at the elevon failure. This is consistent- with the finding through reconligurability calculation [10] that the 
canards are less effective in controlling the pitch movement than elevens, and that loss of eleovn effectiveness 
can significantly affect, tracking the pitch commands. 

The control signals of the LPV and controllers are plotted in Figure 5.3 for the same failure scenario. 
Since the control signals of the LPV controller are the same with or without perturbations, only simulation 




Fig. 5.3. Control signals for each actuator with the LPV and H 0 c . controllers. 

results with “pertl” and “pert2" are plotted in Figure 5.3. It can be observed from Figure 5.3 that the 
elevon signals of the H 0 c controller are significant despite the failure of the elevons, while the elevon signals 
from the LPV controller are insignificant and the canard signals are to compensate the elevon failure in 
Figure 5.3. 

Now, the fault parameters are estimated with an on-line estimator that is integrated into the LPV 
controller as shown in Figure 5.4. The on-line estimator in Figure' 5.4 has two parts: one is a two-stage 
discrete Kalman filter and the other carries out a simple logic of equation (4.7) that converts the bias 
estimates to the corresponding scheduling parameter estimate. The following set of parameter values are 
used in the the two-stage Kalman filter. Sampling time is set at 0.01 se e to capture the response details of 
the open-loop dynamics of the vehicle. The covariance matrices Q J k , Q y and B? k described in Section 3 are 
set as constant matrices with values: 

Ql — 3diag([l, 0.01 2 , 0.01 2 , ().01 2 |), 

Ql — 3diag([Q.Q5 2 , 0.05 2 ]), 

Rl = 3diag([QM 2 , 0.01 2 ]). 

The covariance matrices affect the convergence of the estimator and the noise level considered in the control 
synthesis in Section 4. The initial values of estimated state f 0 |o and biases 7 0 |o are set, as [0 0 0 0] 7 and 
[0 0] T The initial covariance matrices Pj j 0 and PJj 0 are set at IO/2 and IO/4. It is found that the estimates 
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Fig. 5.4. Simulation block diagram of a closed-loop system with an on-line estimator. 

Table 5.1 

Different sets of the varying forgetting factors 


Case 

Ao 

(*min 


1 

0.90 

10 

100 

2 

0.95 

10 

100 

3 

1.00 

10 

100 

4 

1.00 

10 6 

10 7 


are sensitive to the selection to its initial values, but insensitive to the selection of its initial error covariances. 

The most delicate part of the bias estimation lies with the selection of values for Ao, a min and a mQX in 
equation (3.7). This is done by experiments and with little theoretical guidance. Different sets of values in 
Table 5 have been attempted. The fault parameter on-line estimate results are shown in Figure 5.5. It is 
obvious that bias estimation results vary in different cases studied. From the top plots of Figure 5.5, it is 
noticed that the value of X Q affect on convergence rate of estimation. When X 0 is set as 1 at Case 3, the 
bias estimate is not convergent at canard failure situation. For this case, the two fault parameter estimates 
are strongly coupled that can cause false identification of faults. In Cases 1, 2 and 4, initial transients in 
estimates are visible, since the canards are less effective in controlling pitch angle. It is unknown how the 
control surface effectiveness is directly related with the transient behavior of the parameter estimator of the 
two-stage Kalman filter. For Case 4, the covariance matrix P 7 is immediately high value after one step 
integration since a min is defined as 10 6 . It is founded from the results of Case 4 that the high value of the 
covariance matrix leads to good estimate of the scheduling parameter. 

Simulation results at Case 1 are shown in Figure 5.6. “TV” and “VFF’ in Figure 5.6 denote that the 
LPV controller is evaluated at the true values (TV) of the failure parameters and at the estimated failure 
parameters with a variable forgetting factor (VFF), respectively. It is important to note that the LPV 
controller evaluated at the estimated parameter can achieve the desired performance of tracking the pitch 
commands. The difference in tracking performance between using the estimated parameter and the true 
parameter is very small at the steady state. The time delay and transient in the estimate has not formalized 
in the LPV control synthesis process. We currently rely on the robustness of the LPV controller. 

The delay, though undesirable, is helpful in satisfying the rate bounds on the scheduling parameter, 
which is one of the assumptions of the LPV control synthesis methodology. Large delays in fault parameter 
estimates can be detrimental to the stability of a closed-loop system. This is a subject of future study. 
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Fig. 5.5. Scheduling parameter estimates for the cases m Table 5. 


6. Conclusion. In this paper, the LPV controller is designed based on the estimated scheduling param- 
eter. The system variations due to actuator failures are modeled as functions of estimated parameters and 
bounded parameter estimation errors. An LPV controller synthesis problem with bounded parameter errors 
is formulated into two LMI optimizations: one is an LPV control synthesis problem with fixing a scaling 
factor on a uncertainty block and the other is to find a scaling factor given a control law. The optimization 
problem is solved by an iteration method. The performance level of the c losed-loop system with the designed 
LPV controller is reduced bv the iteration method. 

The iteration approach of the LPV synthesis methodology is applied to control of the vehicle at actuator 
failure cases. It is assumed that the actuators are failed one at a time. Actuator failure parameters of 
the vehicle are estimated as biases using an augmented Kalman filter. The LPV controller evaluated at 
the estimated failure parameter is simulated with the plant model which varies as true values of the failure 
parameter. The simulation results show that the LPV controller achieves the desired performance level of 
tracking pitch angle commands for actuator failure cases and robustly stabilizes the vehicle. 
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